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Abstract: The LHC has strongly constrained models of supersymmetry with traditional 
missing energy signatures. We present a variety of models that realize the concept of Stealth 
Supersymmetry, i.e. models with i?-parity in which one or more nearly-supersymmetric parti- 
cles (a "stealth sector" ) lead to collider signatures with only a small amount of missing energy. 
The simplest realization involves low-scale supersymmetry breaking, with an i?-odd particle 
decaying to its superpartner and a soft gravitino. We clarify the stealth mechanism and its 
differences from compressed supersymmetry and explain the requirements for stealth models 
with high-scale supersymmetry breaking, in which the soft invisible particle is not a gravitino. 
We also discuss new and distinctive classes of stealth models that couple through a baryon 
portal or Z' gauge interactions. Finally, we present updated limits on stealth supersymmetry 
in light of current LHC searches. 
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1. Introduction 



Supersymmetry (SUSY) is one of the most compelling possible explanations for the stability 
of the hierarchy between the weak scale and Planck scale, and a leading contender for a new 
principle of nature that could be discovered at the Large Hadron Collider (LHC). The LHC 
is operating successfully, with over 5/fb of data recorded per experiment and many analyses 
of 1/fb of data so far, which have substantially altered the viable parameter space for many 
models of new physics. In the case of SUSY, already in early summer 2011 ATLAS set limits of 
above 1 TeV on gluino and (light generation) squark masses using only 165/pb of data [1]. This 
1 TeV exclusion marked a major milestone on the path to exploring new physics at the weak 
scale, and a plethora of new results followed. Recent theoretical analyses [2, 3, 4, 5, 6, 7] have 
helped to map out the excluded region of parameter space, and have shown that (depending on 
assumptions about which states are degenerate and how they decay) light stops and sbottoms 
may already be excluded up to masses of about 300 GeV. Because the stop cancels the largest 
divergence in the Higgs mass, this indicates that the LHC is poised to begin making definite 
statements about naturalness in standard supersymmetric scenarios. 

As emphasized in [3, 5], two types of studies are currently setting very strong limits on 
supersymmetric models. The first are searches for jets and missing transverse momentum [8, 
9, 10], which suppress standard model (SM) backgrounds by requiring very large missing 
transverse momentum and several very hard jets. These analyses demonstrate the raw power 
the LHC achieves simply by operating at energies never probed by any previous collider. The 
second type of study that sets very strong limits focuses on clean signals that have almost 
no SM background. The exemplar is a search for same-sign dileptons [11, 12], which is a 
very interesting probe of final states with multiple top quarks (among others); for a more 
specialized class of models, diphoton searches [13, 14] play a similar role. The combination of 
jets and missing with same-sign dileptons is already enough to exclude large parts of the 
parameter space of natural SUSY, assuming imparity and decays to invisible particles that 
escape the detector. 

The present stringent LHC limits leave three options for SUSY theories. The first is that 
SUSY is natural, with stops canceling the largest loop corrections to the Higgs mass, and im- 
parity and traditional collider signatures are present. To evade bounds, this requires flavored 
SUSY breaking, with squarks of the first two generations much heavier than those of the 
third. The second option is that SUSY is less natural, and all squarks including the stops are 
heavy, with associated fine-tuning in the Higgs mass. 1 The third case is that SUSY is natural, 
and the squarks of the first and second generation remain light, but the collider signatures 
are altered so that the current searches are evaded: in particular, the classic missing energy 
signature does not apply. One traditional example of the third class is imparity violation 
(RPV). However, as we recently suggested, the third possibility could also be realized in 
SUSY theories with i?-parity, using the mechanism of Stealth SUSY [16]. Stealth SUSY 

1 In some scenarios, a large higgs quartic coupling could be generated, which ameliorates the fine tuning. 
This can be consistent with the possible higgs discovery near mj ~ 125 GeV [15]. 
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opens up the possibility of natural SUSY with light stops and a simple flavor structure. In 
this paper, our goal is to look at a broader class of stealth models. Originally, we considered 
only models with one singlet chiral superfield S as the stealth sector, interacting with the 
MSSM via superpotential couplings to either Higgses or vectorlike matter, and decaying to 
a gravitino. Here, after reviewing the stealth mechanism, we will present several variations. 
We will drop the assumption of low-scale SUSY breaking and explain how high-scale SUSY 
breaking, especially anomaly mediation [17, 18], can also lead to stealth phenomenology. We 
will also consider models in which the stealth sector contains fields with baryon number, 
which lead to events with extremely high jet multiplicities. Furthermore, we will discuss 
some models with new gauge interactions, either a Z' or a confining sector that dynamically 
generates massive singlets. Finally, we will also present an updated discussion of the collider 
physics of stealth SUSY in light of current LHC searches. 

Before ending the introduction, we want to comment on the implication for stealth SUSY 
of the recently observed possible Higgs bump at 125 GeV [19, 20]. Stealth SUSY itself could 
be completely natural with a light stop and small radiative correction to the Higgs mass. To 
raise the Higgs mass above the Z boson mass, additional tree-level contributions to the Higgs 
mass need to be incorporated into the stealth scenario. In this paper, we will not discuss how 
to raise the Higgs mass to 125 GeV, as it could be easily achieved by adding an additional 
module to the stealth SUSY framework. Interestingly, one stealth SUSY model which could 
accommodate a 125 GeV Higgs has already been constructed in [21], realizing the "stealth 
stop" scenario. 

2. General Features of Stealth SUSY 
2.1 The Stealth Mechanism 

The basic ingredients of stealth SUSY are illustrated in Figure 1. The key requirement is 
a set of particles that are nearly degenerate with their superpartners (with supersymmetric 
masses much larger than SUSY-breaking splittings). We will refer to the complete set of 
such fields that feel only small SUSY-breaking as the "stealth sector." It may be as simple 
as a single chiral superfield, as in examples discussed in Ref. [16], or it could be a rich sector 
with one or more gauge groups and many matter fields. In any case, there must be a portal 
through which the lightest (i?-odd) MSSM superpartner (lighest ordinary superpartner or 
LOSP) can decay to a particle in the stealth sector. After this, a decay chain within the 
stealth sector can occur, but it must end with a massive i?-odd stealth particle decaying to 
a nearly degenerate i?-even state plus a light i?-odd state. In the simplest realization, this 
final -R-odd state is the gravitino, but we will be interested in more general models. Finally, 
i?-even stealth states produced in the decay chain must in turn be able to decay back to SM 
fields. The outcome should be that missing energy is carried away only by the light i?-odd 
particle terminating the decay chain, which has momentum suppressed by the small phase 
space available in the decay that produced it. (General decay chains in a complex enough 
stealth sector could involve multiple such lightest R-odd particles escaping the chain; as long 
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Figure 1: A schematic of the sectors involved in a general stealth model. Flavor-blind mediation 
gives rise to standard MSSM soft SUSY-breaking terms, but the soft terms in the stealth sector are 
suppressed relative to this. The MSSM and the stealth sector are weakly coupled, and the size of soft 
terms in the stealth sector is suppressed relative to the supersymmetric mass scale of the stealth sector 
by a weak-coupling factor. 

as the splittings are sufficiently small and the typical multiplicity is low, SUSY can still be 
hidden at colliders.) 

2.2 Stealth SUSY Is Not Compressed SUSY 

It is well-known that, for standard gravity- mediated MSSM spectra, collider signals are more 
difficult to observe as the masses are compressed. For instance, a gluino decaying to a bino 
and two quarks, g — > qqB, is most constrained if the bino is nearly massless, in which case 
a significant fraction of the gluino's energy goes into invisible momentum from the bino. As 
the mass splitting is reduced, the typical missing energy in the event is reduced, and limits 
from LHC searches grow weaker. Recent discussions of limits on compressed scenarios can 
be found in [22]. Superficially, stealth SUSY might sound like a special case of compressed 
SUSY: mass splittings are small, missing Et is reduced, and limits are weaker. However, 
there is a crucial kinematic difference, associated with the fact that in standard compressed 
SUSY, the invisible particle is a heavy decay product, whereas in stealth SUSY the invisible 
particle is very light. This ensures that the reduced missing Et of stealth SUSY is much 
more robust against effects like initial state radiation. 

To clarify this difference, we will review some basic relativistic kinematics and rules-of- 
thumb for hadron collider physics. First, consider the decay of a heavy particle of mass M to 
a particle of mass m = M — 8M and a massless particle. In the rest frame, the momentum 
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of the daughter particles is 



M 2 



— m 



2 



SM 



5M 2 



(2.1) 



V = 



2M 



2M 



Thus, in the limit of a small mass splitting, the momentum of the daughters is approximately 
equal to the mass splitting. Of course, the energy of the heavy daughter is then approximately 
equal to its mass, which is of order M, so the heavy daughter moves slowly, with velocity 
5M/M. The result is that when boosting to the lab frame (with Lorentz factor 7), the 
momentum of the heavy daughter is approximately equal to that of the parent particle, 
whereas the massless daughter will have momentum of order 7 SM, since all of its momentum 
components in the parent rest frame were of order 8M. 

What are the typical (transverse) boosts involved in going to the lab frame for SUSY 
cascade decays? Even without running a simulation or numerically integrating over PDFs, 
we can give a semi-analytic answer to this question. Consider, as an example, production 
of gluino pairs at a proton-proton collider. There are two competing effects that determine 
the typical momentum of a gluino. The first is phase space, which grows with increasing 
momentum; the second are parton luminosities, which determine how often hard enough 
quarks or gluons can be extracted from the protons, and fall rapidly at large momentum. 
Locally, we can approximate that the parton luminosities fall like power laws with exponents 
that will typically be between -3 and -6, depending on details of which partons are colliding 
and in what energy range [23]. In fact, one obtains a reasonably good model of the shape of 
the produced gluino px distribution, da/dpx, simply by multiplying a factor for the growth 
of phase space above threshold (oc px) by a power-law falloff arising from integrating over 
parton distributions (oc (p T + m 2 )~ k ), which amounts to assuming a constant matrix element 
in addition to a simple power-law for the parton luminosities. We give an example in Figure 2. 
A more detailed explanation of why this simple ansatz for the pt distribution is so numerically 
accurate can be found in [24]. The upshot, however, is that the average pt of a produced 
gluon can be found by maximizing the function x(l + x 2 )~ k , where x = px/m; numerically, 
k = 6 works fairly well and leads to a typical px ~ 0.3m. (The distribution is skewed to the 
right, so the mean pt is somewhat, but not dramatically, larger.) 

This means that in the case of compressed SUSY, the invisible particles are not really soft 
at all in the lab frame! For a 600 GeV gluino, for instance, the typical px of a bino arising 
from the gluino decay is 200 GeV, which would be a rather large amount of missing energy. 
As explained clearly in [25, 26], the reason why limits get weaker here is that when gluinos are 
pair-produced, their transverse momenta balance, so that in the compressed case their bino 
daughters also have nearly equal transverse momenta and the net missing Ex is small, even 
though the individual escaping particles are hard. This reduced missing Ex is not robust, as 
in the presence of ISR jets, the gluino pxs are no longer equal, and sizable missing Ex can 
be produced. Thus, compressed SUSY weakens the limits from standard searches, but still 
leads to many events with large missing Ex that can be searched for. 

In the case of stealth SUSY, on the other hand, the invisible particles are soft in the 
lab frame, and the suppression of missing Ex results not from a cancellation but from the 
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Figure 2: Momentum spectra in compressed theories. At left: standard compressed SUSY, with 
nearly degenerate gluino and bino and the decay chain g — > qqB. The bino momentum is typically 
very close to that of the gluino, and is not soft. The orange dotted curve is a simple ansatz da/dpx oc 
Pt(Pt + f" 2 ) 6 to illustrate the characteristic interplay of phase space and steeply-falling parton 
luminosities. At right: stealth SUSY, with the same gluino mass, now decaying in the chain g — > gS, 
S —> SG, and S —> gg. Note that the gravitino, the invisible particle in the stealth case, has a pt 
distribution resembling that of a quark in the usual compressed SUSY case, and is very soft. 

complete absence of high- momentum invisible particles in the event. In particular, because 
the typical transverse boost of the original parent particle (gluino, for instance) is not large, 
we can estimate the boost of the stealth parent (singlino S, in the models of [16]) to be 
7 ~ rrig/mg. Then the lab- frame momentum of the invisible particle is 

Pinvis ~ 7 SM ~ m-g^ 2.2 

m s 

Compared to the bino momentum in the compressed case, which was ~ 0.3 m g , this can 
be made arbitrarily small by taking the stealth splitting small. The reduced missing Ef in 
the stealth case is much more robust, as it is independent of any amount of radiation or the 
structure of the cascade decay. We illustrate some of the relevant px spectra in Figure 2. 

2.3 Stealthy SUSY Breaking 

Having argued that the stealth mechanism is robust from the standpoint of suppressing miss- 
ing energy, the next general issue is whether it is robust from a model-building point of view. 
The setting in which stealthy physics arises with the least effort is low-scale SUSY breaking, 
which always has a light gravitino that appears in the decay of a particle to its superpartner. 
Furthermore, the low scale of SUSY breaking can explain why dangerously large soft terms 
in the stealth sector are absent. One still has to explain the supersymmetric masses in the 
stealth sector, which are near the electroweak scale either by accident or through common un- 
derlying physics. The simplest explanation is to generate them in the same way one generates 
the MSSM //-term; however, to preserve stealthy physics, one would then need to require that 
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Bfi be small. Ordinarily, in gauge mediated model-building, one declares the n/Bfj, problem 
to be solved in models that naturally generate Bfi ~ [i 2 , so we need a more stringent criterion. 
Nonetheless, palatable solutions exist with Bfj, <C fJ- 2 , including retrofitting [27], messenger 
models with i?-symmetries [28], or models similar to the Kim-Nilles mechanism [29] in which 
fi arises from a term of the form (S) n /M n ~ l [30, 31]. In the latter case, one must be slightly 
careful, e.g. arranging for S 2 S 2 /Mp terms in the superpotential to be allowed but not simple 
SS masses; discrete gauge symmetries may be arranged to do this. One could also consider 
stealth sectors with dynamically generated masses (e.g. from confinement) instead of (i terms 
in the superpotential. For low-scale breaking, then, one can formulate robust stealth models 
with the same tools already in use in GMSB models. 

However, one of our main concerns in this paper is to show that the stealth mechanism 
can apply not just for low-scale SUSY breaking, but also for high-scale SUSY breaking. This 
case requires more care, and we will see that the models are necessarily more elaborate. The 
basic model-building problem here is that the SUSY-breaking splittings in the stealth sector 
should be smaller than around 10 GeV. Thus, if m 3 / 2 ~ 10 GeV, we will need to assume that 
generic Mp;-suppressed operators involving stealth fields are absent; that is, we will have 
to assume some form of sequestering [17]. The simplest realization of this is to sequester 
both the MSSM fields and the stealth sector together; perhaps they are localized in an extra 
dimension, far from the source of SUSY breaking. Whether such sequestering is viable in 
string constructions remains a topic of ongoing research [32] , but the problem is no worse for 
stealth models than for more standard phenomenological models, so we will simply assume 
in this paper that it can be achieved. 

Even in sequestered models, there is a remaining technical challenge that will arise in any 
scenario with gravitinos above the 10 GeV scale. It is roughly analogous to the version of the 
n/Bfi problem arising in anomaly mediation models. Namely, suppose that there is a term 
in the superpotential with a dimensionful coefficient, like mX\X 2 . Then we should consider 
m to come with a power of the conformal compensator <p = 1 + 2 m 3 / 2 , such that 

CD J d 2 6 m (1 + 9 2 m 3/2 ) X X X 2 D m 3/2 mX 1 X 2 . (2.3) 

This is a soft SUSY-breaking B-teim given by the mass times m 3 j 2 . Without using the 
conformal compensator formalism, one can also see such terms directly in supergravity. There 
will be terms in the superpotential with expectation value Wq = m 3 / 2 Mp in order to cancel 
the cosmological constant, and cross-terms in the —^2 I W| 2 an d diWdiK-^- parts of the full 
supergravity potential will contain the i?-terms in question. The B-terms are very dangerous 
for stealth SUSY. For example, if we have a mass mS 2 , where S is a stealth field intended to 
have a supersymmetric mass of 100 GeV and a fermion/scalar splitting of 10 GeV, we have: 

5m = m — a/ m 2 — B m — (2-4) 

2m 

Then if we have B = m 3 / 2 m, we require m 3 / 2 < 25m < 20 GeV. This tells us that in any 
scenario in which the gravitino is heavy compared to the stealth SUSY splitting, we must 
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ensure that all of the supersymmetric masses for stealth sector fields do not arise from dimen- 
sionful parameters in the superpotential. They could, for example, arise from dynamically 
determined VEVs, FI terms, or confinement. 

2.4 Portals 

Any stealth model must have some way for the lightest i?-parity odd MSSM particle to decay 
to the stealth sector, and in return for i?-even stealth states to decay back to SM particles. 
A variety of "portals" exist that can allow these decays. Some of the options are: 

• Neutral portals: In this stealth sector operator couples to an operator not 
charged under any continuous symmetry of the MSSM. For instance, given a gauge sin- 
glet chiral superfield S in the stealth sector, one can have in the superpotential \SO neut , 
where O neu t is a gauge singlet MSSM operator. Example portals include O neu t = H u H<i 
(the "Higgs portal"), or O n eut = YY where Y,Y are in the 5,5 representations of 
the SU(5) GUT group. These portals have been discussed in Ref. [16]. Other portals 
include kinetic mixing of U(l) gauge groups, or bifundamental matter charged under 
both an MSSM gauge group and a stealth gauge group. Many more couplings through 
higher-dimension operators are available. 

Neutral portals might be further subdivided; for instance, S could be a singlet under all 
symmetries, in which case it will generically have tadpoles which can be problematic in 
certain models. Alternatively, it might be charged under some nonanomalous discrete 
symmetries, which can help protect against dangerous divergences as well as explain 
the presence or absence of various terms in the superpotential. 

• Charged portals: The stealth states could be charged under symmetries of the MSSM. 
We don't want them to be charged under gauge symmetries, as gauge interactions would 
then typically mediate large enough SUSY breaking to produce missing Et signals. 
Thus, the reasonable candidates are baryon and lepton number. Portals are associated 
with i?-parity-violating terms in the MSSM, with the lowest-dimension examples being 
SLH U , Sudd (the "baryon portal"), SQLd, or SLLE. Such models suggest the pos- 
sibility of unifying stealth collider signals with recent attempts to use such portals to 
explain the abundance of dark matter, baryons, or both [33, 34, 35]. The lepton por- 
tals generically give rise to final states containing neutrinos, potentially reintroducing 
missing Et, so we will limit our discussion in this paper to the baryon portal. 

Notice that in order for the stealth suppression of missing energy to work, we need the 
proper signs of SUSY-breaking splittings in the stealth sector, which depend on which portal 
we wish to use. For example, the portal SH u Hd will always lead to the LOSP decaying to 
one or more SM particles plus the fermion S, so we need rrig > rug- (More generally, the 
scalar and pseudoscalar in S may be split, and we only need one of them to be lighter than 
S.) On the other hand, the operator Sudd lead to the LOSP decaying to SM particles plus 
the scalar S, and in this case we will require that the fermion S is lighter. 
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Aside from a choice of portal, we require a mechanism explaining why a decay of a 
stealth particle to its nearly-degenerate superpartner and a light i?-odd particle will always 
(or almost always) occur in any decay chain. Aside from the gravitino, a Goldstone fermion 
(superpartner of a Goldstone boson) can be a generic candidate for the light i?-odd field. 
More general models will often require a symmetry explanation of why certain decays are 
forbidden; in the case of the baryon portal, baryon number can be useful for this purpose, 
whereas for neutral portals one might consider discrete gauge symmetries. 

Now that we've given a coarse-grained overview of the requirements and challenges in- 
volved in building general stealth SUSY models, let us take a closer look at some particular 
beasts from the zoo of possible models. 



3. Stealth Models with Singlets 

In the section, we will discuss a class of models where the stealth sector consists of singlets, 
either fundamental (in Sec. 3.1) or composite fields (in Sec. 3.5). For the most part, we will 
assume that they communicate with the MSSM through SM vector-like messenger pairs, e.g., 
transforming as 5 + 5 under the SM SU (5). 

3.1 The Vectorlike Portal 

In this section we first review the model in [16] where the singlet S communicates with the 
MSSM through messengers Y, Y, which transform as 5 + 5 under the SM SU(5). 

W D^S 2 + XSYY + M Y YY. (3.1) 

Notice that S must be a pure singlet; any symmetry under which it is charged is broken when 
A and My are both nonzero. Integrating out Y and Y with just the supersymmetric mass 
My at one loop yields the operator 



c J d 2 9SW a W a + h.c, (3.2) 



where W a is the gauge field strength. The matching coefficient c is 

^ lQnMy' v 1 

i 

where the sum runs over the contributions to the gauge coupling beta function coefficients 
from all heavy fields running in the loop. In the paper we adopt the convention that the 
kinetic gauge field term is normalized to — \Tx{F^ u F^ u ) . Writing Eq. 3.2 in components, one 
get operators such as \ a a^ v G aiiV s and sG a ^ u G atJjU . These interactions would induce decays of 
the gluino to singlino plus gluon and of the scalar s to gluons. Similar operators between S 
and other SM vector multiplets exist, which allow decays of neutralinos (charginos) to singlino 
plus 7/Z (W) and of s to two 7's. These decays are illustrated in Figure 3. 
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Figure 3: The decay of a bino LOSP to a singlino and photon through the vectorlike portal, and the 
decay of the singlet scalar to gluons. 



The Y and Y contribution to the strong coupling beta function coefficient is 1. Thus the 
decay widths of the CP even/odd scalar in S to gluons are 



\ 2 a 2 s m\ 
" {a) ~ 64vr 3 m y ' 



(3.4) 



where a s is the strong coupling. Similarly one could obtain the partial widths to two pho- 
tons. Such decays are prompt at colliders if My is near the TeV scale (the lifetimes become 
comparable to those of B mesons or r leptons, i.e. of order 100 /im, when the Y mass is 
about 10 TeV). 



3.2 Tadpole Problems 

For the superpotential we have presented in 3.1, any possible global symmetry under which 
S is charged is broken by the coupling A. Thus S will obtain a tadpole at one-loop order, 
which is logarithmically divergent. Given the superpotential 3.1, at one-loop order, there are 
three tadpole diagrams: two involve a scalar (Y or Y) running in the loop while the third one 
has the ipY,ipY fermions in the loop with a supersymmetric mass insertion. The net result, 
Vtadpole = TS, is logarithmically divergent by dimensional analysis. An explicit calculation 
confirms that 

where A is the messenger scale where SUSY breaking effect is mediated to the low-lying states. 
The tadpole will induce a VEV for S. Consequently, the masses of both fermion and scalar 
components in Y(Y) are shifted: 

m^pY = + rriY 
m Y = (A(5) + my) 2 ± Xm(S) + m 2 soft . (3.6) 

Given A = 0.2, m 3 = 300 GeV, m L = 200 GeV, A = 100 TeV, my = 1 TeV, m =100 GeV, 
we have (S) = 813 GeV and 5m^ Y = 162 GeV. One interesting feature is that, because the 
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supertrace in the Y sector is unchanged, the (S) dependence cancels out of the soft mass for 
the S scalar: 

|A| 2 ,„„ ..ox. Ml 



w 2 ~ 



(6m|) + 4mf,) log 



(3.7) 



(47r) 2 v ' ~ my 

On the other hand, we have omitted terms that are allowed in the superpotential, such as 
kS 3 . In the presence of such terms, a VEV for S could lead to splittings that alter stealth 
phenomenology. Hence, we must assume either that k is small or the portal coupling A is 
small, where "small" in practice means 10~ 2 . 

Thus, in gauge mediation scenarios, the S tadpole is not a severe problem for model- 
building. However, the tadpole can be more troublesome in certain scenarios. Suppose that 
we wish to consider, instead of a decay S — > SG, a decay where another fermion replaces the 
gravitino. One might consider a model 



W 



m 



'-S z + XSYY + M Y YY + -S 2 N. 



(3.8) 



In such a model, in the presence of an S tadpole —TS, the potential is: 



V D \ySN + mS + XYY\ 2 + 



-S 2 



TS. 



The tadpole leads S to get a VEV, and then ./V feels a nontrivial potential: 



dV 
dN 



yS (ySN + mS + XYY 



t 



(N) 



m 

y ' 



(3.9) 



(3.10) 



In this case the VEV for N cancels the leading supersymmetric mass term for S. Other 
models can alter the details, but don't change the basic fact: when S is a pure singlet, it 
inevitably gets a VEV, in the presence of which an S 2 N superpotential term for S decays 
will tend to generate a VEV for TV which spoils the stealth mass spectrum of S. Hence, the 
S 2 N -type models are viable only if S and N can be charged under a symmetry that forbids 
these tadpole problems; we will discuss such a theory in Section 4. 

In fact, the situation is even worse in theories of high-scale SUSY breaking. Naively, 
one might expect that tadpoles are generated only after SUSY-breaking and are at most of 
order m 3 oit . However, as shown in [36], in general tadpoles are quadratically divergent, and 
can be of order m 2 ^ 2 A 2 /Mp in a theory with cutoff A and higher-dimension Kahler potential 
operators suppressed by Mp. Taking A to be of order the Planck scale, we expect that singlets 
will obtain tadpoles of order j^m^^Mp, possibly with an extra loop factor suppression in 
theories with extra structure like no-scale breaking. For low-scale SUSY breaking (gravitinos 
not much heavier than 1 keV) this is not a problem for stealth structure, but for high-scale 
SUSY breaking it is devastating. 

These considerations tell us that, if we want to build stealth singlet models with the same 
phenomenology as that of Section 3.1 that are compatible with high-scale SUSY breaking, 
we should charge S under a symmetry that can forbid it from obtaining a tadpole. 
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3.3 Viable High-Scale Model 

We will now consider an explicit model of high-scale SUSY breaking, in outline. (The full 
details may be found in Appendix A.) As we have just argued, we would like to charge 
the field S under a symmetry to avoid tadpole problems. We will take this to be a discrete 
gauge symmetry. We further wish to avoid large 5-terms, so we would like to start with 
renormalizable superpotential terms and generate masses dynamically. 

The idea, then, is to begin with a superpotential of the form ySQQ + XSYY + kS 3 , where 
the Q fields are charged under a new gauge group and the Y fields are 5 and 5 fields as before. 
The gauge group under which the Q fields are charged should have a number of colors and 
flavors that lead to the dynamical generation of an Affleck-Dine-Seiberg superpotential [37], 
W cS = A fc /det(QQ). Then, as in [31, 38], the singlet fields S and the meson fields Q%Qi will 
dynamically get VEVs related to the scale A. 

In fact, in this model it is possible to arrange for a mesino state, a composite of the 
fields Q, Q, to be the lightest i?-odd fermion, and for a singlino from the S field to decay as 
S — > SM. This decay is prompt on collider timescales. SUSY breaking can be arranged to 
arise from anomaly mediation. The explanation for stealthiness, then, is just that the AMSB 
soft terms are proportional to small couplings among the fields we have introduced. 

The full model-building requires a bit more care in defining an anomaly-free discrete 
symmetry, and a particular hierarchy of couplings, y <C k <C 1. The details have been 
sequestered in Appendix A. 

3.4 Goldstone Fermion Models 

One natural possibility is that the fermion in the final state of the stealth decay is light 
because it is the superpartner of a (pseudo-) Goldstone boson, i.e., it is a Goldstone fermion. 
Perhaps the most familiar examples of Goldstone fermions are axinos. The effective theory 
of Goldstone fermions has recently been discussed in [39, 40, 41, 42]. We consider a (dimen- 
sionless) chiral superfield A with a transformation A—> A + iO. Then in the Kahler potential, 
we can consider terms like: 



where / is the decay constant of the Goldstone. If the broken symmetry is in fact anomalous, 
so that the Goldstone is an axion, then one might have additional couplings that break the 
shift symmetry explicitly: 



Note that both types of couplings, c and b, are of the form that allow a particle to decay 
to its superpartner and a Goldstone fermion. Because we are considering high-scale SUSY 
breaking, we must ask whether it is really natural for the fermionic partner of a Goldstone 
boson to be protected by the same shift symmetry that protects the Goldstone boson itself. 
As recently emphasized in [40] , in the presence of generic Planck-suppressed operators in the 




(3.11) 




(3.12) 



- 12 - 



superpotential one expects the Goldstone fermion mass to be at least of the order of the 
gravitino mass m 3/ / 2 - For high scale SUSY breaking, this is much larger than the stealth 
splitting we are interested in. However, this is not a special property of Goldstone fermions, 
but the generic problem that the stealth sector must be protected from large SUSY breaking 
effects. Thus, we should only consider models in which the Goldstone fermion (and the 
symmetry breaking sector it descends from) are sequestered from SUSY breaking, so that the 
leading dangerous operators in the Kahler potential are absent. 

In order for decays to be prompt on detector timescales, we cannot have / too large. In 
particular, 3.11 corresponds to an interaction 

-jsH A a^d^ s + c.c. = -c^SH^ + c.c, (3.13) 

where the equality holds on-shell in the presence of a mass term msSS. This leads to a decay 
width of the singlino to scalar plus axino which is given in the limit of massless Goldstone 
fermion and small stealth splitting (mg = m$ + 5m) by: 

T { S ^ SA )=^ ms \T) =H WGevUGevJ {— ~ ) 2W ^ 

A QCD axino with / ~ 10 9 GeV and a relatively large coupling c > 1 could be the invisible 
particle A. Or, any symmetry broken at a scale below 10 8 GeV could give a viable candidate 
for this Goldstone fermion. Notice that the phase space suppression in the stealth limit for 
the Goldstone fermion scenario is more mild than in the case of a gravitino, where additional 
powers are present because the interaction is suppressed by the small stealth SUSY breaking. 
One also does not want the decay constant to be much smaller than 10 5 GeV in the presence 
of an AWu coupling, to forbid dangerous decays like g — > gA or B — > ^A that are more rapid 
than decays into the stealth sector. 

3.5 Vector-like Confinement Models 

Now we consider a more complicated hidden sector charged under a nonabelian gauge sym- 
metry. In the setup, there are pairs of vector-like superfields Y and Y transforming as 5 + 5 
under the SM SU(5) and also as 2 + 2 under an additional gauged SU(2)h- The SU{2)h is 
more strongly coupled than the SM gauge groups, and thus we can view the SM as weakly 
gauging the flavor symmetry of the hidden sector. The matter fields Y and Y have a large su- 
persymmetric mass M. Below the scale M, they could be integrated out and we are left with 
a pure SU(2)h super Yang-Mills theory. This model might be thought of as a supersymmetric 
vectorlike confinement model [43]. 

If we assume the 5-plets have a uniform mass M and plug the 1-loop anomalous dimension 
into the NSVZ beta function, we find that at the fixed point, = so the coupling is 
somewhat strong. After integrating out all the fields, the beta function coefficient for pure 
SU(2) SYM is b = 6, and the confinement scale is A = Me" 8 ^/^ 2 ) = Me" 5 / 4 w 0.3M. 
If we want singlet glueballs at the bottom of the SUSY spectrum, near 100 or 200 GeV, 
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this would imply light colored particles from the 5-plets, which decay only through GUT- 
suppressed operators and thus lead to long-lived, R-hadron-like phenomenology. However, 
the CMS collaboration already set a lower limit of 620 GeV at 95% C.L. on the mass of 
a (semi-)stable stop [44]. We expect such a limit also applies here. We could either raise 
the stealth supersymmetric scale a bit higher or split the triplets and doublets in the 5- 
plets (giving them different supersymmetric masses), which could allow further separation 
between the triplet states and the light glueballs. For instance, if the doublet mass Mi is 
300 GeV and the triplet mass M3 is 1 TeV, the one-loop estimate of the confinement scale is 
A = M 3 e" 5/4 (M 2 /M 3 ) 1 / 3 k, 192 GeV. (However, the coupling is already g « 5 at M 2 , so the 
perturbative estimate is likely subject to large uncertainties.) 

In the SUSY limit, an effective action analysis [45] tells us that in the low energy spectrum, 
there are two massive chiral supermultiplets with masses M + and M_ . Both of the masses are 
of order the SU(2)h confinement scale A, which we take to be of order a few hundred GeV. 
According to [45], the heavier one consists of two real scalars, which are SU{2)h gaugino- 
gaugino bound states (gluinoballs) in the limit of no mass mixing term in the potential and 
a Weyl fermion, a gaugino-gluon composite. The lighter one has two real scalars, which in 
the same limit, correspond to parity even/odd glueballs and another fermionic gaugino-gluon 
composite. 1 

SUSY breaking is mediated to the hidden sector through the SM fields. Given an SM 
gluino mass M3, Y and Y receive soft masses at one loop while the SU(2)h gauginos obtain 
soft masses at the two-loop order, 

g 2 glM 3 M mess 

m ^(i6^ l0g — ' (3 - 15) 

where 53 is the SM SU(3) C coupling and M mess is the SUSY breaking messenger scale. In 
the low energy spectrum, this will lead to SUSY-breaking splittings within the gluinoball and 
glueball supermultiplets [47]. Now, in the lighter multiplet, the scalar glueball will be the 
lightest state, with a heavier fermion and a pseudoscalar as the heaviest component. The 
splitting would be of order and is naturally very small. 

Loops of the Y and Y fields serve as portals for the SM to decay to the hidden sector and 
vice versa. For instance, if the gluino is the LOSP, it could decay to the hidden gaugino-gluon 
composites at one-loop. If the gravitino G is the true LSP, the hidden gaugino will decay 
at tree- level to G + gh, or in the composite degrees of freedom, the gaugino-hidden gluon 
composite will decay to the hidden glueball plus gravitino. The hidden glueball could decay 
back to two SM gluons through a box-diagram mediated by Y (Y) fields. 

lr The analysis of [45] is not supported by a recent lattice simulation [46], which orders the glueball and 
gluinoball supermultiplets in the opposite way. The low-energy spectrum of pure super Yang-Mills remains as 
an open question. In this paper, we will stick to the analytic analysis in [45] and [47]. 
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4. Stealth Through The Baryon Portal 



4.1 The Portal: Decays to and from the Stealth Sector 

In this section we will consider the case that a field in the stealth sector carries baryon number, 
and couples to a baryonic operator in the superpotential: 

W D ^ Ui d 3 d k S (4.1) 

Because the udd coupling is antisymmetric in down quark flavors, the portal can have the 
quantum numbers of the A baryon (uds); we will refer to it as the "baryon portal." Note that 
S has charge 1 under U(1)b- As we will see shortly, to have stealth phenomenology the scale 
M is not extremely high, so we will also be interested in UV completions of this operator. 
We will consider either the [/-model: 

W v d a jk Udjd k + MUU + mUuiS + mSS (4.2) 

or the D-model: 

W D D b jk D Uj d k + MDD + b.DdiS + mSS. (4.3) 

Here U and D have the gauge quantum numbers of u and d, respectively, but baryon number 
2/3, and U and D complete them into vectorlike fields. The a's and 6's are flavor-dependent 
coupling constants. Notice that the superfield S contains an R-odd scalar and i?-even fermion. 
In every case, operators containing three quark superfields implicitly have color contracted 
with an e-tensor. We have also added a field S, with baryon number —1, in order to give S 
a supersymmetric mass. For high-scale models, we will want to replace m with a dynamical 
value {X) where the dynamics that give X a VEV may be similar to the model giving rise to 
the S VEV in Section 3.3. 

Our first concern is the condition on A and M necessary for decays to happen within the 
detector. If the LOSP is a squark, it could decay to the S scalar plus two jets: for example, 
one might have £r — > bsS* . Other LOSPs will decay through off-shell squarks to three-jet 
plus scalar S final states, e.g. B,g — > U{djd k S*. We then assume that the scalar S decays 
to its fermionic partner and a light soft fermion (which will be discussed in more detail in 
the following subsections). The fermionic ips then decays back through the portal to three 
jets through a squark-gluino loop. In order to have viable stealth phenomenology, we require 
that both the LOSP decay into S and the tjjg decay back to the SM are prompt enough to 
not create E T signals at colliders. 

To be pessimistic about the lifetimes involved, let us assume that the LOSP is a bino. 
Then it has a 4-body decay B — > Uidjd k S, mediated by a dimension-7 operator: 

£ c s ~ \jk9' ,/ 2 Buidjd k S + h.c. (4.4) 



- 15 - 



Figure 4: Left: one diagram contributing to the decay of a bino LOSP through the baryon portal. 
Right: a diagram contributing to the loop-level decay of the baryon-charged singlino to 3 quarks. 



One of the diagrams contributing to this decay is illustrated in Figure 4. From this we 
estimate 

r ~ g — m 7 - (4 5) 

However, this is very approximate: aside from order-one factors there are functions of rriq/m^ 
and ms/mg that are important. Thus, we have numerically calculated the lifetime, which is 
plotted in Figure. 5. (For this calculation we implemented the interactions, with appropriate 
Lorentz and color structures, in Python in UFO format and used MadGraph5 / MadEvent to 
compute the widths [48].) If the bino is relatively light, say below 300 GeV, and squarks are 
above 1 TeV, the decay is moderate displaced, of order centimeters, for A = 1 and M = 100 
TeV. It can be made much more prompt by considering much heavier binos, but in that limit 
the production of superpartners is likely out of reach of the LHC. Thus, we will consider 
M/A 55 100 TeV as an upper bound on the scale suppressing the interaction, in the case of a 
bino LOSP. 

Less pessimistically, one can consider a squark LOSP decaying directly to the S and two 
quarks. The squark decay width (for one flavor choice) is: 

m - d jdk S) = |gf,/ , (4.6) 



f{x) = 1 + 6 (x 2 + x) log x + 9x - 9x 2 - x 3 . (4.7) 

The phase space factor is important. In particular, if mg = rria/2, we have /(1/4) ~ 0.073, so 
the decay rate quickly becomes very suppressed relative to the case of massless S. Numerically, 
keeping nig = rriu/2, we have: 



(1 Tg\^ "\ ' 
— ) 



M \ 2 1 



100 TeV J A 2 ' ^ 4 ' 8 ^ 



Thus, a squark LOSP decay is prompt. 
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Figure 5: Contours of constant lifetime (in mm) of a bino LOSP decaying through the baryon portal. 
The assumptions are that X u d s — 1 with other couplings turned off, the singlet scalar mass is 100 GeV, 
and AI = 100 TeV. The squark masses are degenerate, ~ rnj = toj, and given on the vertical axis. 



The decay out of the stealth sector is a V>s decay through a squark-gluino loop, illustrated 
on the right-hand side of Figure 4. We find for the decay width: 

r = „ ' ?% L 2 , 4.9) 
3(4^) 5 M 2 m? ' V ' 

with A the dominant and 

1 fl-X 1 



L= dx dy (4.10) 

x + (1 — x)m~/m~ 



o 



assuming squarks are degenerate in mass. The decay length is thus, assuming squarks and 
gluinos have equal masses, 



/100GeV\ 5 / M V/ m- q \2 l 

CT ~ 001cm {^r) (loTw) (iw) 3?' (4U) 

This tells us that, independent of the LOSP, we should not raise M/X far above 100 TeV, and 
even this is slightly delicate, needing a singlino heavier than 100 GeV or gluinos well below 1 
TeV. 

Because the scale suppressing the new physics is 100 TeV or lower, and the couplings 
necessarily couple different generations of down (s) quarks, we should ask whether there are 
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Figure 6: The diagram that contributes to K° — K° mixing in the [/-model. A similar diagram exists 
in the -D-model, along with other diagrams involving Ws or winos. 

strong constraints from flavor physics. The situation is very similar to that of R-parity 
violation with the udd operator, except that we have not actually violated baryon number, so 
processes like nucleon-antinucleon oscillations are not allowed. We expect that K — K mixing 
is the most important constraint. In the ?7-model, the contribution is shown in Figure. 6, 
and is similar to RPV contributions discussed in Refs. [49, 50]. (The D-model has additional 
contributions, but they don't change the qualitative conclusion.) Comparing to the bound on 
the corresponding operator {df^^ s _r)(<^-R7 h s r) as quoted in Ref. [51], we find a constraint: 



This is compatible with our requirements on the lifetime given the above estimates, and 
suggests that M should be near 100 TeV if all couplings are order one. (Note that the A 
couplings are, e.g., X u db oc adbO-u, so the bound doesn't directly constrain them, although 
we should take a u small relative to 4ir.) The bound can be avoided entirely if, for example, 
ads > adb,sb- 

One assumption that might be questioned is the natural size of the couplings \ijk- This 
depends on unknown UV physics. If one assumes MFV, for example, then we should have 
that Xijk oc ^abc{Xu)ai{Xd)bj{Xd)ck up to some unknown constant determined by microscopic 
physics. Assuming the constant is of order 1, even the largest A will be small, of order 10 -3 or 
less [52]. However, since our A's arise at a scale near the TeV scale from products of couplings 
with different flavor structure, it is not clear that MFV is a reasonable hypothesis. Besides, 
even in the MFV framework, the constant could be large, e.g. ~ 10 3 , without violating any 
principle. Detailed considerations of flavor physics and the origin of these couplings is beyond 
the scope of this paper. 

4.2 Baryon Portals with Invisible Gravitino 

Now that we understand how decays into and out of the stealth sector work, we come to the 
problem of the stealth decay itself. The simplest case is to follow the pattern of the models 



M 



Z 160 TeV. 



(4.12) 
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discussed in Ref. [16] and consider the decay S — > ^5^3/2! which occurs promptly for low- 
scale SUSY breaking. For this decay to happen, we require a positive soft mass fh 2 s . In the 
context of low-scale SUSY breaking, one usually considers models of gauge mediation (though 
other possibilities exist). S couples to fields with gauge interactions through either JJuiS or 
DdiS Yukawa couplings, so we can ask whether either of these will automatically produce 
the correct sign for the S soft mass. The usual RG formulas for soft masses generated from 
Yukawa couplings tells us that S will obtain a positive soft mass if + m\ < or, in the 
other model, rh 2 ^ + m 2 , < 0. Within the framework of general gauge mediation (GGM) [53], 
we can obtain different soft masses for U and u only in the presence of an effective FI term 
for hypercharge. In GGM with U, U or D, D and an effective FI term £y , there are two sum 
rules (compared to only one in GGM with £y and without the new fields). One can show 
that in the [/-model, achieving rh 2 ^ + fh 2 ^ < implies: 

rh 2 Q -2rh 2 d -rh 2 L >0, (4.13) 

which cannot be satisfied simultaneously with the usual sum rule 

6m 2 Q -9m 2 d + 3m 2 u -6m 2 L + m 2 e = 0, (4.14) 

without making at least one of d, u, or e tachyonic. On the other hand, in the D-model, the 
condition for rh 2 ^ + rh 2 d < implies: 

ml-rh 2 Q -ml>0, (4.15) 

a condition which can be satisfied. Thus, although we have not specified a full model, the 
GGM formalism tells us that it is reasonable to expect models where the S scalar is heavier 
than the S fermion, and that they will be associated with effective hypercharge FI terms and 
heavy sleptons. (The sum rule shows that when this is true, one also has fh 2 > 9m 2 , + 3m„.) 
This may be an odd corner of parameter space, but there is nothing obviously wrong with it. 
(Indeed, the apparent preference of this GGM analysis for light squarks is nicely consistent 
with the sort of spectrum that provides reasonably prompt decays in our model.) One could 
also consider moving beyond the GGM framework, e.g. by making some SM chiral superfields 
composite as in Refs. [54] and related work. Such models can have different spectra of soft 
masses that are not flavor- degenerate and evade the GGM sum rules. 

Aside from the sign of the S soft mass squared, we also have to check its size. But, as 
we saw in Ref. [16], a splitting appropriate for stealth phenomenology (e.g., ms — m^ s ps 10 
GeV) is easily achieved from RG running (here, between the messenger scale and M) if the 
Yukawa coupling (fe.;, in this case) is ~ 0.1. This is perfectly consistent with the requirements 
for prompt decays and reasonable flavor physics, with M ~ 10 TeV and bjk ~ 1. 

4.3 Baryon Portals with An Invisible Baryon-Charged Fermion 

One interesting model is for S, instead of decaying to its fermionic partner plus another 
weakly-interacting fermion, to decay to a new fermionic field in the stealth sector that also 
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carries baryon number. This leads us to the superpotential (working in the [/-model for 
concreteness) 

W = ajkUdjdt + MUtJ + ai UuiS + aXSS + A5 2 iV 

ra ^ +msS s + xs 2 N (4.16) 

where the U(1)b charge of S 1 is +1, of 5 is —1, and of JV is —2. Notice that the mass term for 
N is forbidden by baryon number symmetry and thus N is naturally light. We have replaced 
a mass term msSS with a Yukawa coupling a to a field X that we assume obtains a VEV. 
This choice is made to avoid the S-term problem of large scalar soft masses ~ ms rn 3/2'i the 
dynamics which gives X a VEV may be similar to that discussed in Section 3.3, relying on 
gauge theories to generate a dynamical scale. We will simply assume that (X) ^ and that 
the X field obtains a mass and does not appear in the decay chains we discuss. Thus, the 
effective theory at low energies is given by the second line of 4.16. 

Because N is a light field carrying baryon number, it potentially leads to rare baryon 
number violating processes in normal matter that could have been detected. One concern is 
double nucleon decay, from n + n — > N' + K° + K° or n + n — > iV' + (7 or 7r°). Diagrams 
contributing to such processes can be found by studying the similar baryon-number- violating 
process in RPV models [49, 55] and attaching S and N lines to the graphs. Two such diagrams 
are shown in Figure 7. 




Figure 7 : Diagrams contributing to double-nucleon decay processes. At left, a contribution to n+n — > 
+ K° + K°, a channel that is open if m N < 884 MeV. At right, a contribution ton + n->JV t +7, 
open if ttin < 1879 McV. Notice that the latter process requires insertions of flavor violation and 
left-right squark mass mixings. Both diagrams require an insertion of a Majorana gaugino mass. 

The n + n ->• iV 1 " + K° + K° decay goes through a dimension- 10 operator, which we 
estimate (up to an order one factor) in the limit in which ^> Mg 3> 7715: 

~ ^r^WM2 a ^ lo s Hr ( uds ) ( uds ) N - ( 4 - 17 ) 

\bTt z m% M?jM z mi, 
s R y b 
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In the case of just an imparity violating \ U( i s coupling, a crude estimate of the bound from 

double nucleon decay is [55] (taking the squarks and gluinos to be at around 1 TeV) ^ 

/ »\5/2 

3 x 10 I 100 MeV/A I , where A is a hadronic scale arising from nuclear matrix elements 



and other order-one slop in the calculation. Relative to that calculation, the decay amplitude 
for us is suppressed by a loop factor and A 2 (and by the fact that this is a 2 — > 2 rather 
than 2—7-3 process), and proportional to the ^4-term and a logarithmic factor. Taking 
A = 10 TeV and a\ = 10 GeV, we can obtain a similar estimate in the stealth case: X u d s & 

5 ^100 MeV/Aj . The hadronic scale A appearing in this estimate need not be precisely 
identical to the one in [55] , but we expect it to be of similar size. If it is 100 MeV, there is 
really no bound at all; if it is somewhat bigger, there may be a very weak bound on A, but 
not one that conflicts with having prompt enough decays for stealth phenomenology unless 
A> 1 GeV. In that case, we would require either that some other flavor choices of Xijt are 
dominant or that the scalar N mass is large enough to forbid the decay. 

In the model as written, the N field is massless. One could add a baryon-number violating 
mass mN 2 , which breaks baryon number to Z/12. This could be taken small, but heavy 
enough that N does not influence BBN, for instance. Such a choice is technically natural if 
this mass is the largest source of baryon-number violation. The remaining Z/12 symmetry is 
enough to forbid proton decay and neutron-antineutron oscillations. 

Now, let us discuss SUSY breaking effects. Again, we need a positive scalar soft mass- 
squared for S for stealth phenomenology, which could be achieved with low-scale SUSY break- 
ing as discussed in Section 4.2. Alternatively, one could try to build a model using anomaly 
mediation, similar to that of Section 3.3. It turns out that the minimal attempt to do so 
yields a spectrum that is only partially stealthy, as discussed in Appendix B; slightly more 
complicated models would be necessary to give viable high-scale SUSY breaking models of 
the baryon portal with S 2 N term. Of course, one can also consider models with alternative 
decays like the Goldstone fermion discussed in Section 3.4. 



5. Stealth Models with a Light Z' 



In this section, we will consider two types of models in which the stealth sector is charged 
under a new gauged U(l)' symmetry with a gauge coupling ga- In the first model, U(l)' is 
broken by an F-term potential. In the second model, U(l)' is broken by a D-term potential 
involving a non-zero FI term. 

5.1 U(l)' Broken by a Superpotential 

5.1.1 Decays to and from the Stealth Sector 

In addition to the MSSM superpotential, our model has the following superpotential, 



W = yS(<M_ - A 2 ), 



(5.1) 



where the chiral superfields have opposite charges q, —q (q > 0) under U(l)' while S 

is not charged. The scale A will be taken to be around the weak scale and could be generated 
dynamically through the retrofitting mechanism [27]. 

For this sector to communicate to the MSSM, (some) SM fields (and possibly SM-charged 
spectators to cancel anomalies) are also charged under this U(l)'. One would then imme- 
diately worry whether such a low scale Z' sector is already ruled out experimentally. First, 
to avoid stringent constraints on leptonic processes, this Z' has to be leptophobic or have 
very tiny couplings to the leptons. For a leptophobic Z' , the main constraints come from 
the dijet resonance searches at hadron colliders and electroweak precision tests (EWPT). The 
strongest limits on a Z' decaying into dijets with mass below 500 GeV actually come from a 
combination of UA2 [56] and CDF results [57] (LHC dijet searches are less sensitive to low 
mass resonances). At the Tevatron's energy </s = 1.96 TeV, a Z' produced at a rate of a 
few tens of picobarns is allowed across the whole range below 500 GeV and a rate of hun- 
dreds of picobarns is allowed below mz> = 350 GeV. For an order one coupling ~ 1, the 
charge of the first generation quarks under this U(l)' has to be small < 0.1. This is naturally 
achievable if the first generation quarks only obtain a non-zero charge of order ~ 0.1 after 
mixing with a heavy vector-like set of quarks with order one charge under the U{\)' [58]. 
This effective Z' point of view also relieves us of the need to carefully construct anomaly-free 
symmetries. Aside from direct constraints, Z 1 gauge boson mixing with the SM Z boson 
provides a positive contribution to the oblique T parameter and thus is constrained by the 
EWPT. Parametrizing the off-diagonal component in the 2x2 mass matrix of Z and Z' as 

2 

7m|, the correction to the T parameter is then aT = 7 2 2 f 2 with a the fine structure 



constant. Fixing = 125 GeV, mt = 173 GeV and U = 0, the global electroweak fit gives 



S = 0.07 ± 0.09 and T = 0.10 ± 0.08 [59], which implies 7 , ™ z < 0.036. If our Higgs 



is charged under this U(l)' , 7 = 9d9 ^ v where g is the SM SU(2)w coupling and EWPT 
requires the Higgs charge under U(l)' Qh < 0.03 for mz' = 200 GeV. 1 In summary, the 
experimental data only constrains the charges of the first generation quarks and the Higgs. 

At the scale A, U(l)' is spontaneously broken by the VEVs of the <j) fields = (4>-) = 
f. Ignoring SUSY breaking for the moment, / = A. All chiral superfields (S, </>+, and 
(/>_) obtain masses. The fermion in the linear combination 4>\ = (<fi + + 0_)/v2 obtains a 
Dirac mass mjy = y2yf with the fermion in S, from the superpotential. The orthogonal 
combination <p2 = (<A+ — <P-)/V^ is "eaten" by the super-Higgs mechanism to make the whole 
gauge multiplet massive with mass m g = 2gaqf ■ More specifically, the U{1)' gaugino Z' 
marries the combination (f>2 = (4>+ — <^_)/v2 to obtain a Dirac mass m g . The imaginary part 
of 4>2 is eaten by the gauge boson while the real part of (p2 gets a mass m g from the D-term 
potential. 

A squark LOSP charged under U(l)' could decay to the quark and the U(l)' gaugino, 
q — > qZ' , while a [/(l)'-charged Higgsino LOSP could decay to the Higgs and the U(l)' 

1 For a two Higgs doublet model with a large tan/3, the down type Higgs could have a charge of order one 
as its contribution to the mass mixing of Z and Z' is suppressed by tan 2 f3 [60] . 
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Figure 8: Possible stealth decays through the Z'-portal. The decay on the left passes through the 
Z', while the decay on the right passes through the scalar within the massive vector supermultiplct. 

gaugino, H — > hZ . Other LOSPs not charged under the U(l)' could decay through charged 
MSSM particles, e.g. off-shell squarks, to two jets plus a Z . After SUSY breaking which we 
discuss in the next section, a small mass splitting is generated between different components 
of the heavy vector multiplet. Depending on the split mass spectrum, Z' could decay as 
Z 1 — > Z' + X with X the soft invisible fermion or to Re(c/>2) via the mixing with the (fo, 
Z' — > Re(02) + X. The Z' would subsequently decay back to two jets, Z' — > qq. For the 
scalar Re(02), the trilinear scalar coupling from the D-tevm would induce its decay through 
an off-shell squark-gluino loop to two quarks, Re(</>2) — > qq. These possible decay paths are 
illustrated in figure 8. All decays in and out of the stealth sector are prompt. 

5.1.2 SUSY Breaking: Decays inside the Stealth Sector 

Now we consider the effects of SUSY breaking in the stealth sector. We will proceed in three 
steps: we will first parametrize SUSY breaking by three phenomenological parameters, the Z 
gaugino soft mass m^, and the ) soft masses and show that they have to satisfy 

certain conditions to fulfill the stealth mechanism. Then we will check how to satisfy these 
conditions in different SUSY mediation mechanisms. Finally we discuss the decay inside the 
stealth sector during which a light invisible particle carries away missing energy. 

After SUSY is broken, the whole scalar potential is V = \F\ 2 + \D\ 2 + m 2 + \4)+\ 2 + m 2 _\4>_\ 2 . 
For simplicity, we will assume m + = m_ = m s , which is true in a broad class of SUSY breaking 
mediation schemes that are insensitive to the sign of the charge. The minimum is shifted to 
/ = A 2 — m 2 /y 2 . The uneaten scalars' masses are 



m im{4> x ) = V2yA, m s = m Re ^ = y / 2(y 2 A 2 - in 2 ) w V2yA - 




= ^A>-(^-i)*J - W - (f - JL) f , (5.2) 



where we expand to the leading order in fh 2 / A 2 . The gaugino soft mass ra^, would push one 
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fermion mass eigenvalue down and the other one up. Analytically the fermion masses are 

fh?' 

m 4>i = m s = V 2 (y 2A2 - m i) ~ ^2y A - ^/2y A > 

m^ 1(2) »2 gd qA-^-f±-f-, (5.3) 



where ^1(2) are mixtures of Z 1 and fa. The gauge boson mass Z' is 



mz> = 2 9d qf = 2g d qJtf - ™± « 2g d qk - ^. (5.4) 

V y y A 

Neglecting the 5, </>i sector, we see that the mass of the gauge sector, to the leading order in 
rhg/A 2 , is shifted by an overall — after SUSY is broken. Another feature is that the 

Z' mass is always in between the two fermions fa and fa with \mz' — m Vi(2)l ~ l m z'l- ^hc 
lightest state is a boson, Re(fa), only when the following conditions are satisfied: 

I ~ 2 1 

m, < 0, J — 5j- > m-, . (5.5) 

The first condition requires that the scalar get a negative soft mass while the second one tells 
us that the mass splitting between the scalar and fermions is much larger than that between 
the gauge boson and fermions. 

Now we will look for mediation mechanisms that generate a soft mass spectrum satisfying 
Eq. 5.5. One irreducible SUSY breaking contribution to the stealth sector is U(l)' gauge 
mediation through the visible sector. The gaugino Z' receives soft masses at one loop with 
the MSSM fields (or the spectators that cancel the U(l)' anomaly) running in the loop. As 
the gaugino soft mass breaks R symmetry, there must be an insertion of a chirality-flipping 
scalar soft mass, which breaks R symmetry, in the loop. Assuming contributions only from 
the MSSM fields, the Z 1 soft mass is 

3g%q?m t/ n m- g%q?usin2B m\ 

m Z>= 022 ( Va * Sin & ~ ^ COS $ lQ g ~4 ~ 1 a 2 lQ g » 5 - 6 

where v is the Higgs VEV, qt (qh) is the charge of the top quark (Higgs) under the U (1)', and f3 
is the mixing angle between up- and down-type Higgs. fi is the dimension one parameter in the 
MSSM Higgs superpotential and ttia is the mass of the heavy Higgs scalars and pseudoscalar 
in the decoupling limit, at is the top A term and is the stop mass where we neglect the 
mass differences between left and right handed stops. The contributions from other quark 
supermultiplets are negligible as they always proportional to the fermion masses. 

Now we consider the scalar soft masses from Z' mediation. The fields </>+ and <j>- obtain 
equal soft masses at the two loop order, with MSSM fields as messengers. 2 The MSSM 

2 We have assumed here that the one-loop contributions to the scalar soft masses proportional to the trace 
over the U(l)' charged fields, Tr(ij m^), vanish. This is true in general gauge mediation without a non-zero 
FI term. 
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messenger mass matrix has a positive non-zero supertrace, which leads to a negative and 
logarithmically divergent two-loop contribution to the soft masses of both <f> + and 0_ . The 
leading logarithmic contribution from one multiplet labeled by i is [61] 



32vr 4 
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STr m { log — z 



uv 



m- 



(5.7) 



where Ni counts the degrees of freedom for species i, e.g. N = 3 for a particle transforming 
as (3,1) under the SM SU(3) C x SU{2)y/] STrm 2 is the supertrace of the messenger mass 
matrix; m is the mass of the scalar messenger; and the UV cutoff A is the scale where the 
soft masses of MSSM fields are generated. For g& = 1, q = 1, iVj = 3, qi = 0.8, m = 400 GeV 
and Ajjv =1000 TeV, we find m 2 ~ —(40 GeV) 2 . For the supersymmetric mass m g = 100 
GeV, this leads to a splitting of 8 GeV between scalar and fermions inside the stealth sector, 
neglecting m^,. 

According to Eq. 5.5, the gaugino soft mass has to be more suppressed compared to the 
scalar soft mass. In gauge mediation, this is easily achievable as the A terms are generically 
small and m^, calculated from Eq. 5.6 is suppressed. There could be other SUSY breaking 
contributions to the stealth sector, for instance, gravity mediation in scenarios with a heavy 
gravitino. To maintain the desirable stealth spectrum, besides sequestering the scalar soft 
masses m 2 , we require that ra 2 remain negative and the mediation mechanism has to be 
approximately R symmetric. We will not go further in building complicated high-scale me- 
diation models to realize this but only mention that there is no no-go theorem against these 
requirements and there are already viable high-scale mediation models with an approximate 
R symmetry in the literature [62]. 



220 




Figure 9: One sample spectrum of the Z' model. Dashed lines correspond to two fermion states, 
which are mixures of U(l)' gaugino and (0+ — </>_)/\/2 fermion. We only show states in the gauge 
sector, which are relevant for the decays. States such as S and <f> + + do not mix with the other 
states and will not appear in the decay chain. 



Finally we put in the invisible particle and estimate the decay lengths of the decay inside 
the stealth sector. In gauge mediation, the natural candidate for the invisible particle is 
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the gravitino. The decay ^1(2) ~~ > B.e{4>2) + G has a decay length ranging from mm to cm 
depending on the SUSY breaking scale y/~F as long as the splitting between the fermion 
and scalar is about 10 GeV. The heavier fermion ip2 could also decay to the gauge boson, 
V> 2 -> Z'(Z'*) + G. As discussed in [16], in decays A -> B + G, if T(J3) ~ 0(1) GeV, the width 
of the three-body decay through an off-shell B will be comparable to that of the two-body 
decay through an on-shell B particle. As the Z' width could be big, e.g. T(Z') ~ 1 GeV given 
that the charges of all the light quarks are 0.1, the width of ^2 — * Z'* + G is comparable 
to the two-body one and thus this channel is not completely stealthy! On the other hand, 
as the scalar Re(4>2) decays through a loop to two SM quarks and its width is much smaller 
than 1 GeV, the decay tpi/2) Re{(p2) + G is dominantly two-body. Luckily as the mass 
splitting between the gauge boson and gaugino determined by m?, always has to be smaller 
than that between Re(4>2) and fermions, e.g. by a factor around 10, the partial width of the 
channel through gauge boson is suppressed and thus negligible. One natural replacement for 
the light gravitino is the Goldstone fermion, as discussed in Sec. 3.4, as long as it is light 
(say, m£ < 10 GeV) and cosmologically safe. A could couple to the (ft+ and 0_ fields in the 
Kahler potential as K D (A+ A^)<fr i(f>+ + (A-\- A^)(jr_(j)— , and the estimate of the decay length 
is the same as in Eq. 3.14. 

5.2 Light Z' from Kinetic Mixing 

In Section 3.3, we discussed the use of a strongly coupled gauge theory to generate a dynam- 
ical scale. An alternative mechanism for dynamically producing a mass term is to generate 
a VEV using kinetic mixing of U(l) gauge groups, — \e f d 2 9W a W a . This is the super- 
symmetrization of the phenomenon studied in [63]. As shown in [64] (also see [65]), it can 
generate a supersymmetric mass scale in a sector with no explicit superpotential mass terms, 
starting from a U(l) in a different sector with a D-term expectation value. This will produce 
a VEV v 2 = e(.D). We can first ask whether a hypercharge D-term VEV can suffice. The 
experimental limits on kinetic mixing with a new light Z' are typically at least as strong as 
e 55 0.03 for most Z' masses below or of order 100 GeV [66]. Suppose we take, for example, 
Dy = (2 TeV) 2 , e = 0.03, a dark U(l) coupling g& = 0.3, and a dark sector superpoten- 
tial W = \S(j)+(/>- with A = 0.2. Then the D-term potential g\ ^|<^+| 2 — |<^-| 2 — £,dj, with 
£d = eDy, will lead to a VEV Vd = (<f>+), producing multiplets of dark photon states at mass 
9d v d ~ 100 GeV and of dark Higgs states at mass Xv^ ~ 70 GeV. As in [64], we expect a 
scalar dark Higgs to be lighter than its fermionic partner if kinetic mixing is the dominant 
mediation of SUSY breaking to the dark states. (One could also consider AMSB, arranging 
for gd large enough relative to A to produce a similar spectrum.) The dark Higgs decays 
through two off-shell dark photons; using the results in [67], we can see that for the relatively 
heavy dark Higgs mass and kinetic mixings that we consider, this decay is prompt, unlike in 
previous phenomenology. The cascade decays proceed much as in the study of lepton jets [68], 
with a bino decaying to dark Higgs and dark Higgsino, B — > h^h^, followed by — > h^G 
given low-scale SUSY breaking or perhaps to a Goldstone fermion in the case of high-scale 
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Figure 10: In this section, we will consider the LHC limits on the two stealth diagrams shown in 
this figure. The left diagram is an example of fully hadronic stealth SUSY, where the LOSP decays 
to jets and soft gravitinos. The right diagram is an example where one photon is produced in each 
SUSY cascade. 

breaking. Dark Higgses then produce jets of charged quarks and leptons; because of the mass 
scales we consider, these will no longer be primarily lepton jets. Similar decays through a 
dark sector already have been known to hide SUSY; here, the challenge is made even greater 
by combining this with the stealth mechanism. 

Alternatively, one could consider models in which the portal to the MSSM is independent 
of kinetic mixing, but kinetic mixing with a new heavy gauge boson generates a supersym- 
metric mass scale in the stealth sector. For instance, in models of low-scale SUSY breaking, 
it may be that new U(l) gauge groups exist at the messenger scale. Then a relatively small 
kinetic mixing, such as the 1CP 3 to 1CP 4 mixings that we expect to be typically generated 
by loops, can lead to weak-scale FI terms for a stealth sector U(l), inducing VEVs that can 
be used to generate stealth-sector masses. This provides an interesting alternative to the 
nonabelian gauge dynamics needed to generated masses along the lines of Section 3.3, while 
preserving the portals and collider signatures of such models. 

6. Stealth at the LHC 

In this section we discuss several aspects of stealth LHC phenomenology. As we saw above, 
the stealth SUSY framework can lead to many different final states, depending on the portal 
that connects the MSSM to the stealth sector, and on the identity of the LOSP. The common 
features of stealth phenomenology are low missing energy and high multiplicity final states, 
because the LOSP decays through the stealth sector to SM particles. We begin by evaluating 
the limits that the LHC can now set on stealth SUSY. In section 6.1, we discuss the limit on 
gluinos that decay to a fully hadronic final state, as in the left of figure 10. In section 6.2, 
we evaluate the limit on stealth decays that produce photons and jets as in the right of 
figure 10. Finally, in section 6.3, we show that displaced vertices can lead to extra missing 
energy, although the size of the effect is typically small. 
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6.1 Limit on Stealth with Jets and MET 



Here we present the limit on fully hadronic stealth decays. We focus on gluino pair production, 
where the gluinos decay through the process shown to the right of figure 10. Each gluino 
decays to three jets and a gravitino. This is the dominant final state in the SYY model 
with a gluino LOSP. The momentum of the gravitino, and therefore the amount of missing 
energy, depends on the mass splitting, 5m between S and S. We consider the limit coming 
from jets plus missing energy, which will allow us to evaluate how small 5m must be in order 
to sufficiently suppress missing energy and hide SUSY. 



search 


luminosity fb 1 


6m limit [GeV] 


ref. 


CMS jets + MET 


36 x 1(T 3 


17 


[69] 


CMS jets + MET 


1.1 


12 


[9] 


ATLAS 6-8 jets 


1.34 


13 


[70] 


ATLAS 2-4 jets 


1.04 


35 


[8] 


CMS a T 


1.14 


38 


[71] 


CMS M T2 


1.1 


45 


[72] 



Table 1: The LHC searches for jets and missing energy included in our analysis. The "6m limit" 
column lists the value of the splitting in the stealth sector below which the search sets no limit at all, 
for a 100 GeV singlino. We have included all of the jets plus missing energy searches conducted with 
1 fb -1 , except for the CMS search based on the razor variable [73], where we expect similar results. 
For comparison we have also included the 2010 search that we found to set the strongest limit on 
stealth SUSY [16]. 

In our previous paper [16] , we considered the limit on the same process coming from the 

2010 data, L ~ 35 pb - . We found that the strongest limit in 2010 comes from the high 
Ht channel of the CMS jets plus missing energy search [69]. With this search, we found 
no limit when 5m < 17 GeV, fixing the singlino mass to 100 GeV and restricting to gluinos 
heavier than 200 GeV. Here, we consider the limit that is set by various 2011 jets plus missing 
energy searches, all with about 1 fb , listed in table 1, that were carried out by ATLAS and 
CMS. We use pythia 6.4 [74] for event generation, PGS 4 [75] as a crude detector simulator, 
and K factors for gluino production from Prospino [76]. The limit that we find is shown in 
figure 11, as a function of the gluino mass and 5m, for a singlino mass of 100 GeV. We find 
no limit for 5m < 12 GeV, and this comprises the stealth regime where SUSY is hidden from 
searches that demand missing energy. In the large 5m regime, the gluino mass limit extends 
to about 800 GeV. In 2011, as in 2010, we find that the strongest limit on stealth comes from 
the CMS high Ht search strategy, which is not surprising because this search demands less 
missing energy than the other searches, compensating this with a tighter cut on Ht- 

One surprising feature of the limits shown in figure 11 is that, in the small 5m limit, the 

2011 searches present only a modest improvement over the 2010 limit, despite the factor of 
30 increase in luminosity. The reason for this is that the 2011 searches have tighter cuts and 
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LHC Limit 




Figure 11: The limits from searches for jets plus missing energy on stealth SUSY with fully hadronic 
decays, shown as a function of the gluino mass and the mass splitting in the stealth multiplet, 5m. We 
find that hadronic stealth SUSY is completely unconstrained by these searches when 5m < 10 GeV. 
Interestingly, the limit coming from 2011 searches is comparable to the limit coming from 2010 searches 
in the small dm regime. 



therefore significantly lower acceptance than the 2010 searches. The acceptance for stealth 
SUSY to pass these cuts is shown in figure 12, where we compare the acceptance of the CMS 
high Ht searches from 2010 and 2011. The 2010 version of the search demanded missing 
energy above 150 GeV and Ht > 500 GeV. The 2011 version increased the missing energy 
cut to 200 GeV and the H T cut to 800 GeV. 

6.2 Limit on Stealth with Photons 

Now we consider the LHC limit on stealth decays that produce photons. We focus on squark 
pair production, where each squark decays to 1 photon and 3 jets, as to the right of figure 10. 
Every event contains two hard photons whose presence can be used to greatly reduce the 
QCD background. We find that a limit can already be set on this scenario using the 77 mass 
distribution, which has been used by CMS [77] and ATLAS [78] to set limits on KK gravitons 
with about 2 fb _1 . Although stealth SUSY does not produce a narrow 77 resonance, colored 
production can lead to a large enough cross-section to exceed the measured spectrum. 

In order to check the LHC limit, we consider a spectrum where all light flavor squarks are 
degenerate and decay according to the diagram of figure 10. We fix the bino mass to 300 GeV, 
the singlino mass to 100 GeV, and the stealth scalar mass to 95 GeV. The right side of 
figure 13 shows the m 77 spectrum measured by CMS compared to the stealth spectrum that 
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Figure 12: The acceptance for stealth SUSY events to pass the high Ht selection of the 2010 and 
2011 CMS searches for jets and missing energy. These are the two searches that dominate the 2010 and 
2011 limits, respectively. Despite the large increase in luminosity in 2011, we find competitive limits 
in the small 8m regime because the increase in data is compensated by the decrease in acceptance 
that results from the tighter 2011 cuts. 

would be produced by rriq = 500 and 700 GeV. We use Pythia 6.4 [74] for event generation 
and we simulate the photon isolation ourselves from the hadronized Pythia output. We use 
Prospino [76] for the NLO squark production cross-section. To the right of figure 13 we show 
the 95% C.L. limit that we derive on the degenerate squark mass, rriq > 700 GeV. In order 
to derive this limit, we have binned the spectrum in 100 GeV bins, and we require that no 
bin exceed the 95% C.L. limit using the CL S statistic [79]. 

To summarize, we have found that stealth decays to photons are not too stealthy, and 
rriq > 700 GeV is required for degenerate light flavor squarks. However, this limit is weaker 
than the limit on degenerate squarks in regular SUSY from jets plus missing energy, rriq > 
900 GeV, after 1 fb -1 . We believe that there is still reach for a dedicated analysis to discover 
stealth decays to photons because the background can be further reduced by demanding a 
large of jets in addition to the presence of two hard photons. Furthermore, it may be 
possible to reconstruct the 'yjj resonance, as we discussed in our previous paper [16]. 

We also note that certain stealth models include narrow, low mass, 77 resonances, and the 
branching ratio into these resonances can be used to constrain the total SUSY cross-section. 
For example, the vectorlike portal of section 3.1 leads to a small branching ratio of the stealth 
singlet to 77. In the benchmark model of this type from our previous paper [16], the singlet 
branching ratio to photons was 4 x 10~ 3 . A branching ratio of this size is constrained by 
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the 5 fb _1 higgs searches [80, 81] if the singlet mass is within the search region of m s ~ 
100 — 150 GeV. The limit on a x Br is on the order of 0.1 pb, which constrains the SUSY 
cross-section to be smaller than about 50 pb. This requires the gluino to be heavier than 
about 300 GeV. 
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Figure 13: Stealth SUSY decays that produce photons, Yi — > jjj, are constrained by the 77 mass 
distribution. The left plot shows m 77 measured by CMS [77] with 2.2 fb _1 and the mass distribution 
that would be produced by degenerate light flavor squarks that decay to a neutralino that decays to 
ijj- We have fixed the neutralino mass to 300 GeV, the singlino mass to 100 GeV, and the stealth 
multiplct splitting to 10 GeV. The right plot shows the ratio of the signal cross-section to the 95% C.L. 
limit that can be derived from the 77 spectrum. We find that nig < 700 GeV is excluded. 



6.3 Missing at Long Lifetimes 

Here we wish to present an effect mentioned in [16] in more detail. If missing transverse energy 
is computed from energies in the calorimeter, the pseudorapidity and azimuthal angle (</>, 77) 
may be used as detector coordinates relative to the primary vertex to compute the momentum 
of the particle that deposited the energy, assuming it originated at the collision point. This 
will continue to give the correct momentum for particles originating at displaced vertices, as 
long as the decay products continue along the line of flight of the decaying particle. More 
generally, however, it will lead to a mismeasurement of momentum, as shown schematically 
in figure 14. 

To illustrate this effect, we have simulated 300 GeV gluinos decaying to 100 GeV singlinos 
in the decay chain g — > gS, S — > SG, S — > gg, shown to the left of figure 10. We compute 
the distribution of decay locations for the S — > SG step, for various singlino lifetimes cr. 
Assuming a calorimeter located 1.25 meters from the beamline, we then compute the measured 
missing Et based on the locations where gluons from the S decay will hit the calorimeter. 1 
If the scalar S is not extremely boosted, the S — >■ gg decay has a wide opening angle, and the 

A more complete treatment, which we leave for further work, would also need to include the calorimeter 
shower shape, which can change the jet 4-vector for very large displacements. 
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Figure 14: Jets that originate from a displaced vertex will be measured to have the incorrect mo- 
mentum because the momentum of a jet is normally determined by pointing the calorimeter hits back 
to the primary vertex. This momentum scrambling can lead to missing energy in stealth SUSY events 
with displaced vertices. 




Figure 15: The missing Et distribution for different singlino lifetimes, shown separately on linear 
and logarithmic scales. We find that displacement does enhance the missing energy tail although the 
effect is mild, even for displacements of 10s of cm. Here we have fixed the singlino mass to 100 GeV 
and the stealth multiplet splitting to 10 GeV. 



location of the calorimeter energy can vary significantly as the lifetime is varied. In Figure 15, 
we show the resulting missing Et distribution for various cr values on both linear and log 
scales, in the case of a splitting m^ — mg = 5 GeV. The distributions peak in about the same 
place at all lifetimes, but the high missing Et tails become much larger at lifetimes of order 
the detector radius. (Only those events for which the decay happens before the calorimeter 
radius are plotted.) For lifetimes of 10 cm, though, the distribution has hardly changed at all. 
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Figure 16: The difference of the measured missing energy and the Monte Carlo truth missing energy, 
for difference singlino lifetimes. We see that larger lifetimes lead to a larger difference, and that the 
missing energy is more likely to be increased than decreased. 

This is reassuring and in accord with our geometric intuition. A similar illustration appears 
in Figure 16, this time illustrating the difference between the measured missing Et and the 
true value. Since the overall missing energy distribution is not too sensitive to this effect, the 
limit we found in section 6.1 carries over to the case where the singlino has moderate lifetime. 

Of course, if the tracker is used instead of just the calorimeter, events with 10 centimeter 
lifetimes will look dramatically different from those with only prompt decays. We encourage 
searches for macroscopically displaced vertices or anomalous jets without tracks pointing at 
the calorimeter deposits, with all possible trigger paths. 

7. Discussion 

We are now truly in the LHC era, with 5/fb of data accumulated by each experiment and 
numerous beyond the SM analyses completed using 1 to 2/fb. Already, the landscape of 
possible theories has changed dramatically. Gluinos, in most realizations of i?-parity preserv- 
ing SUSY, are excluded up to masses of 700 GeV or more. Squarks of the first and second 
generation are similarly constrained, pointing to either flavored mediation of SUSY breaking 
or tuning in the MSSM. The possible Higgs discovery at close to 125 GeV also points either 
to extensions of the MSSM or to fine-tuning in SUSY. We must be prepared to move beyond 
the MSSM if we are to continue to take SUSY seriously as a natural theory of the electroweak 
scale. 

In light of the emerging LHC data, it is becoming increasingly important to consider 
extensions of the MSSM that are less constrained by the LHC. By surveying such extensions, 
new collider searches can be motivated in order to fill the gaps in existing search strategies 
and lead the way to unexpected discoveries. We have presented a rich class of models that 
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realize this idea through the stealth SUSY mechanism. Missing energy is removed, and in 
its place are high multiplicity final states with rich resonance structures and the possibility 
of observably displaced vertices. If the end of the 7 (or possible 8) TeV run at the LHC fails 
to bring a SUSY discovery in the missing energy tails, then it will become essential to assess 
the experimental status of stealth SUSY. There remain many promising and unconstrained 
channels where natural SUSY may reveal itself. 
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A. Details of the High-Scale SYY Model 

We will consider an explicit model that can work with high-scale SUSY breaking, as briefly 
outlined in Section 3.3. We have argued in Section 3.2 that, to avoid tadpoles that spoil the 
stealth mechanism, we wish to charge our field S under a symmetry. To keep the same spirit 
as the models of [16], we will take this to be a discrete gauge (nonanomalous) symmetry. 
Furthermore, as we argued in Section 2.3, a simple mass term of the form mS 2 or mSS in 
conjunction with high-scale breaking leads to a l?-term soft mass that is not compatible with 
stealth phenomenology. Thus, we would like to write down a model with only dimension-3 
operators in the superpotential, which nonetheless dynamically generates a mass for S and 
spontaneously breaks the discrete gauge symmetry at a scale low enough that the induced 
tadpole is safely small. 

Suppose we introduce a new SU(4) gauge group with 3 flavors Qi, as well as 3 new 
vectorlike fields in the 5 + 5 of the SM SU(5), Yi, Yi. We also introduce three singlets Si. We 
consider a superpotential: 



This superpotential respects an anomaly- free Z/3 discrete symmetry under which all fields 
have charge 1 (see for instance [82] for a guide to discrete anomalies). The choice of W 
breaks the flavor symmetry acting on the Q's and Y's down to a product of U(l)s. This is 
an admittedly ugly tactic to avoid having a moduli space of supersymmetric vacua, because 
the massless modes on the moduli space can be problematic for us. This potential will give 




(A.l) 
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S a VEV and hence a dynamical mass, as in [31, 38]. The SU(4) gauge group confines and 
generates an ADS superpotential [37], 

w <* = d4W ■ (A2) 

In the presence of this dynamical mass scale A, the theory has a supersymmetric vacuum 

(S i ) = (£) 1 \ (A.3) 



(QiQi) = (f) 19 ^- ^ 

In the low-energy theory, instead of working in terms of the Q, Q fields, we will work in terms 
of the meson fields Mij. In particular, we will assume that a canonically normalized meson 
is related to the quark fields by a matching condition: 

QiQj = vkMij, (A.5) 

with v an order-one number. Let us compute the mass spectrum; it is easiest to work with the 
fermions. W e ff gives rise to Dirac masses ^ 2 y 5 ' 3 K _5 / 9 A^My V'-Mj, for the off-diagonal mesinos, 
while the diagonal M fields mix into two states with that mass and one state (Tr M) with four 
times that mass. The singlinos get masses — K 5//9 y 3 / 9 A from the S 3 terms, and also mix with 
the diagonal mesinos through the yvSiMa terms. Our goal for this sector is to generate a 
supersymmetric mass for the singlet fields S, as well as for the Y fields, without substantially 
affecting SUSY breaking for the S fields. 

To calculate SUSY-breaking effects, we assume that the MSSM fields together with the 
S, Y, and Q fields are all sequestered, and that AMSB is the dominant contribution to the 
soft masses [17, 18]. The SU(4) confining sector poses a complication for computing soft 
masses, as do the Y fields, as nonsupersymmetric thresholds can take AMSB off its usual RG 
trajectory [83, 84, 85]. If the new states are below the gravitino mass, one can evaluate the soft 
terms above their mass scale and then run down, taking the thresholds into account. If they 
are heavier than the gravitino mass, they will contribute non-AMSB scalar mass-squareds of 

order 167r 2j^-2 [84]. We will take A 3> m 3/ / 2 so that this is the case for the composite states. 
There are two interesting hierarchies of couplings to consider in which the Si and states 
are not highly mixed. If k < y < 1, the mesinos are heavy relative to the singlinos, and 
{S) S> (M). Because the YY mass is X(S), this implies very small values of A if the Y 
fields are to be light enough for their interactions to contribute any sizable AMSB soft terms. 
As a result, even the suppressed m^ 2 /A 2 contributions from the composite sector can be 
important, and the SUSY-breaking spectrum is not fully calculable. 

The limit y <C k <C 1 is the more viable one, and leads to mesinos much lighter than 
the singlinos. Let us make some numerical estimates in this case. We fix 772.3/2 = 30 TeV 
and A = 125 TeV, taking y = 10~ 5 , k = 10 -3 . (We set v = 1, for the unknown matching 
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of canonically normalized mesons with fundamental quarks; other values could be accom- 
modated by changing other parameters to absorb the difference.) Then the supersymmetric 
mass of the singlet states is 116 GeV, two of the meson states have supersymmetric masses of 
40 MeV, and the other meson state has a supersymmetric mass of 120 MeV. The mixings are 
small, e M s & 10" 2 . The singlet VEV is (S) « 58 TeV. We take the Y fields to be below the 
gravitino mass, my = A (S) = 3 TeV, by choosing A = 0.05. Because A > n,y, the dominant 
contribution to the AMSB soft term for the singlets comes from A and is given by: 

™ 2 s = \ K/ 2 | 2 f t l! = |^ (35 |A| 4 - f |A| 2 , 2 ) - -2400 GeV 2 . (A.6) 

Ignoring all other SUSY-breaking effects, this would correspond to a stealth splitting m§ — 
mg ~ 10 GeV, precisely as desired. (Strictly speaking, this calculated soft mass is valid above 
the scale my, though we expect running below my to make little difference.) There are a 
number of other sources of SUSY-breaking effects to check, however. Given the term —^kS 3 
in the superpotential, AMSB gives rise to a scalar trilinear (A-term) ^tu^^^kS 3 , which in the 
presence of an S VEV becomes a contribution to the scalar mass. However, f3 K oc k\ 2 / (16-7T 2 ), 
and this contribution is negligible. 

The next concern is the contribution from the meson states, which mix with the singlets. 
The mesons feel SUSY breaking in two ways: first, there are AMSB contributions from 
the superpotential term W e g = z/ -3 A 6 / det M; second, there are power-suppressed effects 
from the nonsupersymmetric thresholds at the scale A associated with composite states that 
interact strongly with the mesons. Let's first tackle the AMSB contributions. Expanding 
around the VEV, M = (M) + SM, we have schematically A 6 /detM ~ (A/ (M)) 6 5M 3 . 
Numerically, for our parameters (M) ~ 22A and so the ADS superpotential, for the purpose 
of estimating AMSB soft terms, amounts to a Yukawa coupling of order 10 -8 and makes no 
significant contribution. More significantly, we expect soft masses of order fh 2 M ~ 167r 3 /^2 from 
interactions with the composite states. These should lift the scalar meson states to masses 
~ 500 GeV. However, we expect that the S scalars are insulated from these effects, as they 
do not interact directly with the composite states. Through mixing to a meson and back, we 

estimate a contribution 5fh 2 s ~ e 2 MS 167r 3 /^; ~ (0.3 GeV) ms, negligible relative to the AMSB 
soft term S acquired through its interactions with Y, Y. 

Finally, we observe that this model has a built-in stealth decay of a singlino to its scalar 
partner and a mesino, S — > SM, arising from the S 3 vertex and the small mixing €ms- The 
non-stealthy decay S — > MM is forbidden by additional mixing factors (smaller than the 
phase space suppression) and by the expectation that the scalar M is lifted to larger masses 
by its interaction with the heavy composite states. The stealth decay S — > SM has a lifetime 
of order 10 microns, and would appear prompt from the point of view of collider physics. 

This gives a proof of principle that a stealth model compatible with high-scale SUSY 
breaking is possible. Sequestering and a discrete gauge symmetry enforce the basic structure 
of small SUSY-breaking splittings and absence of dangerous i?-terms. However, we had to 
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choose Yukawa couplings of order and 10 -5 : not smaller than those we have seen in 

nature, but perhaps uncomfortably small. Also, due to confined, strongly interacting states 
at relatively low masses that affect the SUSY-breaking mass terms for mesons mixed with the 
stealth fields S, the model is not as calculable as we might prefer. Finally, the breaking of a 
discrete symmetry raises the specter of domain wall problems in cosmology, much as in the 
NMSSM, and a somewhat more sophisticated version of the model could be needed to evade 
them. Stealth SUSY, then, seems to find a slightly less congenial home in the context of 
high-scale SUSY breaking than it did for low-scale SUSY breaking. It would be interesting to 
explore further whether alternative models might work without such small Yukawa couplings 
and light composites. 

B. SUSY Breaking for the Baryon Portal with S 2 N 

We will now give details of AMSB for the case of the baryon portal, Sudd, along with an 
S 2 N operator to allow decays to a baryon-charged invisible fermion. We will sequester both 
the MSSM sector and the stealth sector. In the low-energy theory below the U or D mass, 
the Yukawa coupling XS 2 N contributes a positive soft mass. However, this is not the full 
answer. The fields U or D may be lighter than the gravitino mass, in which case they give 
nonsupersymmetric thresholds that take AMSB off its usual RG trajectory [83, 84, 85]. In 
that case, we should still apply the AMSB formulas, but only above the scale M, and then 
run the results down from the nonsupersymmetric threshold according to the usual RGEs. 
Above the scale M, new Yukawas also contribute. Thus, we should evaluate the AMSB soft 
masses using all interactions in the first line of 4.16. The result is: 

|2 

a ^20|A| 4 + 12|A| 2 |a| 2 + 4|A| 2 |a| 2 ) (B.l) 



~2 m 3/2' 

m N - 



(16tt 2 ) 



2 



ms = (40 |A| 4 + 36 |A| 2 \a\ 2 + 15 \a\ 2 - 16 \a\ 2 (g 2 + 9 — 



+12|A| 2 H 2 + 6|a| 2 H 2 + 2|a| 4 + 16vr 2 \a\ 2 j$) (B.2) 
2 

4 \a\ 2 |A| 2 + 3 \a\ 2 \a\ 2 + 2 |a| 4 + 16vr 2 \a\ 2 j$) (B.3) 



-2 m 3/2 
m 



S (ia^2\2 



(16tt 2 ) 

am ' i/2 ^4|A| 2 + 3|a| 2 + 2|a| 2 + 167r 2 7 |) (B.4) 



16vr 2 

where |a| 2 implicitly means ^ |aj| 2 and the anomalous dimension 7^ cannot be calculated 
without a UV completion of the mechanism by which (X) is determined. Because X has a 
VEV, the trilinear a a term becomes a B-term: 

BSS = - m ^™2 2 (4 |A| 2 + 3 \a\ 2 + 2 \a\ 2 + 16^ 2 7 f ) SS. (B.5) 

The — \a\ 2 term in m| threatens to give a sizeable wrong-sign soft mass if we do not take 
a <C A. Furthermore, the I?-term splits the eigenstates and tends to produce a scalar lighter 
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than the fermion. This is difficult to avoid, as m| S> \B\ typically implies a large enough A 
that splittings are not stealthy. This differs from the vectorlike portal with AMSB discussed 
in Section 3.3 because in that case we could take the piece of a soft mass proportional to — g% 
to dominate. Here, that would give the wrong sign. 

Anomaly mediation can thus easily give stealth-size splittings, by taking the couplings 
to be small, but it tends to produce well-mixed scalar states, one of which is lighter than its 
fermionic partner. The outcome is a partially stealthy spectrum: about half the time, the 
LOSP will decay to the heavier scalar S, which will decay to its partner S and a soft fermion 
N. The other half of the time, the decay will proceed to the lighter scalar S. Its decay is 
not stealthy, because it will go through an off-shell S (in a four-body decay). Comparing 
to Equation 4.11, we can estimate that for A, X u d s ~ 0.1, ~ 200 GeV, and gluinos and 
squarks at 500 GeV, the lifetime would be of order a few centimeters. In these decays the N 
is sharing energy with several visible particles, and for small splittings the singlino propagator 
is enhanced when the N has minimal energy, competing with phase space suppression in that 
region. Thus, we expect that missing energy is strongly reduced even in the half of the decays 
that are not, strictly speaking, stealthy. 

Alternatively, we can simply appeal to an extension of AMSB to produce soft terms com- 
parable in size that lift both scalar states above the fermion. Because minimal AMSB in the 
MSSM predicts tachyonic sleptons, a number of models have been proposed to generate addi- 
tional contributions to SUSY-breaking for scalars of the same order as AMSB contributions 
without changing the basic paradigm of sequestering (e.g. [83, 84, 86]). Such models can be 
adapted to produce a fully-stealthy sequestered version of the Sudd model with S 2 N decay. 
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